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In the context of superfluid Fermi gases, the Ginzburg - Landau (GL) formalism for the macroscopic wave function has been 
successfully extended to the whole temperature range where the superfluid state exists. After reviewing the formalism, we first 
investigate the temperature-dependent correction to the standard GL expansion (which is valid close to Tc). Deviations from the 
standard GL formalism are particularly important for the kinetic energy contribution to the GL energy functional, which in turn 
influences the healing length of the macroscopic wave function. We apply the formalism to variationally describe vortices in a 
strong-coupling Fermi gas in the BEC-BCS crossover regime, in a two-band system. The healing lengths, derived as variational 
parameters in the vortex wave function, are shown to exhibit hidden criticality well below Tc. 


1. Introduction 


The Ginzburg - Landau (GL) approach is a powerful tool 
for the description of superconductors in the close vicinity of 
the critical temperature Tc. Recently, the GL method was 
re-derived in the context of superfluid ultracold Fermi gases 
iSSlill. The GL approach was also applied to explain the 
phenomenon of the “L5-type” superconductivity. However, the 
validity of the GL approximation far below Tc is still under dis¬ 
cussion 10]. In this connection, much efforts were undertaken 
to extend the GL mproach to a wide range of temperatures (see, 

e. g., RefsJ0ISIlE[iIi[il)- 

In Ref. il3n . we formulated an extension of the GL theory for 
a two-band superfluid fermion system solvable for the whole 
range 0 < T < Tc assuming slow variation of the order param¬ 
eter in time and space, without any assumption on the magni¬ 
tude of the order parameter. The theory is mainly focused to the 
strong-coupling ultracold atomic Fermi gases in the BCS-BEC 
crossover. In the present work, we briefly review the method 
and the description of vortices in a two-band system at temper¬ 
atures T < Tc, where the standard GL technique is apparently 
inapplicable. 


The formalism developed in Ref. Ill311 is aimed mainly at the 
investigation of localized deviations of the order parameters T'j 
from a uniform equilibrium background A,. These deviations 
can be, for example, vortices or solitons il4 1 151 llhl llTl IlSl 
HH. A frequently used theoretical method to study these 
localized deviations at temperatures far below Tc is a Bogoli- 
ubov - deGennes (BdG) equation set. Re-formulations of the 
BdG method for ultracold atoms can be found, e. g., in Refs. 
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1 211 l22l |23L |24|] . The present method can be used as a com¬ 
plementary tool to the BdG equations and is straightforward to 
implement numerically. Moreover, the BdG equations are re¬ 
stricted to the mean-field approach, while the present method 
can be used beyond the mean-field approximation accounting 
for fluctuations about the saddle point. 


2. Formalism 


The starting point of our treatment is the partition function of 
a two-band fermion system in the path-integral representation. 


Z cx 




( 1 ) 


with the fermion action depending on the Grassmann fields 
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-I-/ 


dr U (r, t) . 


( 2 ) 


Here, 5 o is the free-fermion action functional, 

72 


5o = 




H-o-.j I (3) 


accounting for both spin and band imbalance of the fermion 
system - through unequal masses mj and chemical potentials 
p/. The fermion-fermion interaction U is given by: 


+ ^3 (<AT,l'AT,l*Ai,2<A'i,2 + ^i,I<Ai,l^T,2*/'T,2) 

+ ^4 (iAt,1'At, 1^T.2'A'T,2 + ^i,l<Ai,l^i,2|Ai,2) ■ (4) 
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It contains the terms describing both the intraband s-wave scat¬ 
tering (with j - 1,2) and the scattering between fermions in dif¬ 
ferent bands (with j = 3,4). Introducing the auxiliary bosonic 
fields and performing the Hubbard - Stratonovich transforma¬ 
tion we arrive at an effective bosonic action of the pair fields 
as described in Ref. Q. Subsequently, we make the standard 
approximation for the GL approach; we assume that the pair 
fields slowly vary in time and space. The gradient expansion of 
the pair fields leads to the long-wavelength approximation for 
the effective bosonic action. This yields the following GL-like 
free energy 
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In order to analytically compare the results of the present 
approach with the known GL method near Tc, we use the results 
of Ref. which represents the limiting case of the present 
approach when T —> Tc (for a one-band system and without 
imbalance). For temperatures near Tc, the order parameter is 
small. Thus the coordinate-dependent thermodynamic potential 
Q.sj is expanded in powers of |'1',| up to the quartic order, and 
the coefficients Cj, Dj and £, are kept for T', = 0. In this case, 
the GL-like free energy © is reduced to the TDGL free energy 
of Ref. ll^ . except for the coefficient D, which appears to 
be real in the present approach. The reason for this difference 
consists in the following. The imaginary part in D appears in 
Ref. when the gradient expansion is performed at T* = 0. 
On the contrary we perform the summations of the whole series 
in powers of the order parameter before taking the limit T —> 
Tc, indicating therefore that the appearance of an imaginary part 
of O depends on a sequence of the limits 'T —> 0 and T —> Tc- 

3. Results 


Here, the function j formally coincides with the saddle-point 
thermodynamic potential for the imbalanced Fermi gas, 
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where Fkj' = Bogoliubov excitation energy, 

^k.j - Sr ~ free-fermion energy, and the chemical po¬ 
tentials for the imbalanced fermions are /ij - + /r, /2 and 

jl. However, the order parameter 'T, entering 
this thermodynamic potential is coordinate-dependent. The pa¬ 
rameter y describes the strength of coupling between two bands. 
The coefficients C,, Dj and £, derived in Ref. iQ are 
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The functions fp (p, s, f) are the Matsubara sums; 
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They can be analytically expressed, e. g., using the recurrence 
relations; 
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First, we illustrate a difference of the temperature behavior 
of the coefficients of the GL-like free energy © compared to 
the coefficients of the TDGL equation of Ref. 12-511 . 



Figure 1: The coefficient C calculated, as a function of temperature, within the 
extended TDGL formalism (solid curves) and within the TDGL theory of Ref. 
flSl (dashed curves). 


In Fig. [T] the coefficient C is plotted as a function of tem¬ 
perature for several values of the inverse scattering length l/a^ 
and compared with the coefficient c of Ref. 1251] . Both co¬ 
efficients analytically tend to the same values when T —> Tc- 
The temperature behavior of the coefficient C drastically dif¬ 
fers from that for the corresponding coefficient c within the GL 
approach ll^ . In the whole range of the BCS-BEC crossover, 
the coefficient C only slightly varies when T goes from Tc to 
zero. When increasing the inverse scattering length, the range 
of temperatures, where C and c are rather close to each other, 
gradually broadens. In the molecular (BEC) regime, both so¬ 
lutions tend to one and the same limit for all T < Tc- On the 
contrary, for 1/aj < 0, (i, e., at the BCS side and at unitarity) 
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c rapidly increase when decreasing temperature (except for the 
BEC case), and even diverges at T —> 0. These results confirm 
the fact that the standard GL approach becomes inapplicable at 
low temperatures. 

Vortices are studied in the present work using the variational 
method. The deviations of the order parameters 'Ey from a uni¬ 
form equilibrium background can be represented through the 
product of a uniform background amplitude Aj = with 

the amplitude modulation function fj (r, t) and the phase factor 

giSj{r,T). 

= A,-/,(r,T)eW. (13) 

The coefficients Dj,Cj,&j are kept with the bulk values of the 
order parameter. Thus the time and space dependence are taken 
in leading order through the derivatives. This is in line with 
the gradient-expansion approximation which was already kept 
when deriving (|5]l. 

Further on, we introduce the notations; 


j nij 

„„ _ (Cj + 2£)|aJ 


(14) 

(15) 


The parameter is the superfluid density, and is the 
quantum pressure coefficient. 

Using we arrive at the following variational GL-like 
free energy functional. 



I |2 , 90i 



^jm\m2y 
7.71 


A 1 A 2 / 1 / 2 COS (02 - 0l) 


(16) 


It describes, in principle, not only the stationary states but also 
the time-dependent Josephson physics for a two-band system 
due to the phase difference 02 (r, t) - 0i (r, t). 

In the present work, the amplitude modulation function for a 
vortex was used in the form 


/ (r, ^) = tanh 



(17) 


with the healing length The healing lengths are determined 
minimizing the free energy (fThl l. In Fig. |2]a, the healing lengths 
for a vortex in a two-band system are plotted as a function of 
temperature for the inverse scattering lengths (in units of the 
Fermi wave vector kp) l/ai = 0, l/a 2 = -0.5, and for dif¬ 
ferent values of the coupling parameter j. The healing length 
for the “stronger” band, fi, extremely weakly depends on y. 
The healing length for the “weaker” band, ^ 2 , demonstrates 
the “hidden criticality” discussed in our manuscript and in Ref. 
iIt^ . At zero interband coupling, each of two subsystems (the 



Figure 2: {d) Healing lengths for a two-band superfluid fermion system as a 
function of temperature for different values of the coupling parameter y; {b) the 
ratio of the healing lengths ^ 2/^1 • 


“stronger” and “weaker” bands) is characterized by its own crit¬ 
ical temperature Tcj and healing length which tends to in¬ 
finity at T ^ Tcj. When the Josephson interband coupling is 
nonzero but sufficiently weak, we can see a fingerprint of the 
phase transition for a “weaker” band as a peak of the healing 
length ^2 at T » Tc, 2 . 

When comparing the healing lengths for a vortex calculated 
in the present work note with those calculated in Ref. Q us¬ 
ing the model fermion system near a hard wall, we see a qual¬ 
itative agreement between the healing lengths determined by 
these two methods. However, there is some quantitative differ¬ 
ence between these healing lengths. 

In Fig. |2]l7, the ratio ^ 2/^1 is plotted for the same parameters 
as in Fig. 1. We can note on a remarkable similarity between 
these results and those shown in Fig. 2 (c) of Ref. ll^ for a 
two-band superconductor using BdG equations. The ratio ^ 2/^1 
starts from a value ^ 2/^1 > 1 at zero temperature, exhibits a 
peak near the critical temperature for a “weaker” band 7'c,2, and 
tends to 1 when T ^ Tc (which is very close to Tc,\). 


4. Conclusions 


In summary, we re-formulated the path-integral approach for 
interacting Fermi gases ll25n to the case of a two-band system. 
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The Hubbard-Stxatonovich transformation and the integration 
over the fermion fields lead to an effective bosonic action with 
Josephson interband coupling. The gradient expansion of the 
effective bosonic action results in the GL-like free energy func¬ 
tional in which the amplitude of the pair field is not a small 
parameter. Therefore the obtained free energy represents an 
extension of the Ginzburg-Landau formalism to temperatures 
below Tc- The range of applicability of the gradient expansion 
is determined by the same conditions as for the standard GL 
approach, where that expansion is also used. Thus the present 
extended GL-like method is valid under the same conditions as 
the GL approach - but in a wider temperature range. 

As an example, the method has been tested for vortices in a 
two-band system of ultracold fermions. It has been shown that 
the “hidden criticality” far below Tc, treated previously using 
the BdG equations 12611 is captured by the extended GL-like ap¬ 
proach. Because of the validity of the present approach at tem¬ 
peratures far below Tc, it can find a wide spectrum of applica¬ 
tions, e. g., for the analysis of distributions of trapped fermionic 
atoms, vortices, solitons and other spatially non-uniform phe¬ 
nomena in ultracold Fermi gases. 
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